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ABSTRACT 

We examine the 0(a) forward-backward asymmetries for the production process e + e~ — > 
7, Z — > qq(g), tagging the outgoing heavy-quark jet at center-of-momentum energies off 
the Z-peak. The complicated analytic results are reduced to simple polynomial forms 
that provide excellent approximations. For charm and bottom quark, a full dynamical 
cancellation gives 0(a) zeros in the forward-backward asymmetry close to the Z-peak. 
We conclude with a detailed numerical analysis of our results. 



PACS number(s): 11.38.Bx, 11.80.Fv, 14.65.-q 



*Feodor-Lynen Fellow 



The measurement of asymmetries in the production of fermion pairs at e + e~ colliders has 
been proven as an indispensable tool to examine rigorously the most important properties 
of the Standard Model. In general, these experimental techniques have advanced to such 
an extent that theoretical predictions beyond the Born approximation in the perturbative 
series of the relevant couplings have to be taken into account in order to agree with the 
given measurement precision. 

Of particular interest is the forward-backward asymmetry A FB which results from 
the vector (V) /axial- vector (A) interference terms of the intermediate 7, Z bosons in the 
production process e + e~ — ► 7, Z — > //. Therefore, its measurement allows for a direct 
determination of the relative strength between the V- and ^4-components of the fermion 
coupling to the neutral current, or, equivalently, for a precise determination of the effective 
electroweak mixing angle sin 2 ^. Experiments that measure A FB for various quark flavors 
and leptons are carried out by the different LEP and SLD collaborations. 

This work concentrates on dynamical QCD one-loop effects in the forward-backward 
asymmetries of heavy quarks. For the charm and bottom quark, high-precision measure- 
ments on the Z-peak are performed by the ALEPH, DELPHI, and OPAL collaborations 
at LEP In the theoretical literature, a first analytical treatment of massive 0(a) 

corrections for A q FB at the Z-peak can be found in an article by Djouadi et al. M. Only 
recently their result (given as an expansion in the quark mass) has been slightly corrected 
by Stav and Olsen ||. The exact analytical formulas are lengthy and complicated, see 
e.g. Ref. §. 

However, in the following we shall present compact Schwinger representations for 
the C-odd structure function in the differential production cross section for heavy quarks. 
To the best of our knowledge no such representations have been treated in the literature 
before. Subsequently, we use these results to find simple polynomial expressions for the 
single-jet forward-backward asymmetry including 0(a) radiative corrections. These ap- 
proximate formulas give very accurate estimates valid over the entire physically relevant 
energy spectrum, which will then allow us to reveal interesting dynamical properties of 
A q FB at QCD one-loop level. Finally, we conclude this work with a detailed numerical 
analysis for charm-, bottom-, and top-quark production. 
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For heavy-quark production e + e~ — > 7, Z — > qq, the differential cross section is 
usually integrated over the azimuthal angle to yield the following decomposition in terms 
of the polar angle 9 of the scattered quark 

f°" - = i(l + cos 2 0)o- t/ + f sm 2 9a L + % cos 9a F . (1) 
a cos 9 s 

The structure functions a v and a L correspond to unpolarized and longitudinally polarized 
gauge bosons, respectively, and their sum gives the total cross section a T = o v + o L . 
Clearly, in Eq. (U) the term containing a F is the only component that changes sign under 
the replacement 9 — > n — 9, and thus constitutes the odd term under charge conjugation 
C in the fermionic final state. 

At the Born level, a F is straightforwardly given by 

2 

a F/0 = 8n^v 2 g VA , (2) 
q A 

where q is the momentum transfer carried by the exchanged 7 or Z boson, and v = 



1 — Am? /q 2 with quark mass m. The factor g VA incorporates all couplings that result 
from the mixed VA interference in the intermediate state 

g VA = -Q q a e a q ReXz + 2 v e a e v q a q \x z \ 2 ■ (3) 

Here, the fractional charge of the quark is Q q and the relevant electroweak coulings are 
Vf = 2 T/ — 4Qf sin 2 9 W and a/ = 2 T/ for / = e, q. The Z-propagator displays the typical 
resonance behavior for the decay of an instable massive particle so that we have 

XM 2 )= 9 :f lq l * with g F = ——1=- ~ 4.299 ■ 10 5 GeV 2 . (4) 



q 2 - Ml + iM z T z ^ 8V2 



not 



The calculation of o F including first-order corrections in the strong coupling involves 
the summation of virtual soft-gluon contributions and hard-gluon bremsstrahlung. A 
closed form expression was derived in Ref. 0: 



°>/l — °~F/0 



all ~ 
1 + Re A + Re C) v + -± C F - \ - (4 - 5^)S 2 - ^(1 - 0(§3 + §e 



47T V 

2(4 - 30§ 4 + e§6 + 2(§ 8 + S 9 + 3§ 10 - §n) + 2(1 - 0(2 - £)Sia 



(5) 
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where C F = 4/3 is the usual Casimir operator of the SU (3) color group and an additional 



mass parameter £ 



Am? /q 2 has been introduced. 



The coefficients A and C are the conventional QCD form factors in the notation of 
Refs. (using the Feynman gauge and on-shell renormalization), and correspond to the 
gluonic corrections of the V and A currents, respectively. For the axial form factor C, it 
has explicitly been shown that different methods like anticommuting 75 within dimensional 
regularization M, dimensional reduction ||, and the 't Hooft-Veltman prescription for 
75 || all produce identical results. 

In the following, we shall only be interested in the high-energy (v — > 0) and low- 
energy (y — > 1) limits of the form factors A and C 



a s 



1 Re A ~ Re C ~ -^C F 
Air 



ln 2 £- (l + 41n2) Inf + 2(1 + 21n2) In 2 



+|tt 2 -4 + 0(0 







Re A 



ReC 



a i 
Air' 



Air' 



C P 



ir~ 



v 



- 8 + ir 2 v + |(1 + 241n2> 2 + 0() 



(6) 



(7) 



4 + ir 2 v - |(11 - 241n2> 2 + 0() 



Eq. (0) clearly states that for nearly massless quarks or sufficiently high center-of- 
momentum energies, E cms = y/q 2 , the virtual-gluon corrections become insensitive to the 
parity property of the relevant vector-boson vertex. However, in the asymptotic energy 
range near threshold, A and B differ by a finite contribution reflecting the distinct nature 
of the underlying symmetries. Following the reasoning by Schwinger [nj, both 1/f-poles 
in Eqs. (^) and (H) correspond to the strong attraction between the color charges (quarks) 
in the non-relativistic limit with a relative velocity 2v <C 1. 

In Eq. (g) the collinear IR divergences emerge as logarithmic singularities for £ — >• 0, 
and eventually cancel when the hard-gluon parts are added. (Note that the trivial soft IR 
divergences have already been neglected as indicated by the wiggle, viz. Ref. 0.) The full 
analytic solutions of the qqg phase-space integrals Sj, i = 2, . . . , 12, have been calculated 
and classified in Ref. 0. Much simpler results which approximate the exact solutions in 
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the important high- and low- energy domains were presented in Ref. [O] [[ 

It is now straightforward to take the limits in the total expression Eq. (||). For 



1, we obtain 



a,0( 



(9) 



i. e. the QCD corrections for a F are genuine quark-mass effects and can safely be neglected 
in the high-energy region. This prediction was already made in Ref. |§. On the other 
hand, the asymptotic behavior close to threshold has never before been considered in the 
literature. Our result for v — > is 

^2 1 



1 + ^C F --6 + A + 0Q 
2tt v 



(10) 



Although non-perturbative resonance effects at threshold supersede predictions made by 
perturbation theory, Eq. flT0D gives essential information on o F sufficiently above the 
production resonance. 

Combining Eqs. (|9|) and ( |I0"D yields the following Schwinger representation for o F 

^ = l + C F a s U 1 -- V (v)), (11) 
a F /o 2 V v ) 

with the mass-zero condition <p(l) = 1. This mass-zero condition is an absolute require- 
ment for any ^-representation, whereas the exact threshold value y?(0) = 6/ it 2 is of lesser 
importance due to the 1/f-pole dominance in Eq. (Pp. In general, <p(v) is a function of 
< v < 1 and contains apart from the universal term in Eq. (|TID (describing the color 
interaction of the quarks close to threshold) all the non-trivial energy dependence of <j f/1 . 
Suitable (^-representations that provide excellent approximations to the exact solutions 
are simple polynomials of degree m > 2: 



in 



E CLiV 1 

<Pm(v) = ^ • (12) 

E a-i 

i=0 

Tab. 1 displays the coefficients for the lowest-order representations <p2, (f3, and (p^. 
In Fig. 1, we have plotted these polynomial representations together with the exact re- 
sult. Already the second-order form <£> 2 gives a very accurate interpolation so that the 
t Note the typographical error in Table I of Ref. jy]. The correct limiting behavior of the 'spin- 
dependent' integral 85 close to threshold (v — > 0) is §5 ~ 4( 21nw — 1 ) + 0(). 
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corresponding Schwinger formula Eq. ([□]) provides a very compact expression for easy 
implemention of the 0(a) corrections to a F . 

A straightforward procedure to include these massive 0(a) effects already in the 
Born approximation consists in the introduction of effective couplings. From Eqs. (||) and 
(PI) we obtain directly the prescription 

7T / 1 



9 



VA 



~9 VA = g VA 



1 + C F a s - l - - r \„{r) 



(13) 



where (p m (v) are the appropriate polynomials of Tab. 1. 

Similarly, one finds replacements rules for the two remaining couplings, g and 
g AA , that multiply with the C-even components of the differential rate. They contribute 
through the VV and AA parity-parity combinations of the intermediate 7, Z states to the 
to the total rate 



(7 7 



Anct 
q 2 



1 /o 2\ VV 1 3 

±v (3 — v ) g + v g 



3„AA 



(14) 



Using the third-order Schwinger representions given in Ref. ||11|| , we find the following 
explicit 0(a) expressions 



~vv 
9 



x 



Q q -2Q q v e v q ReXz + (v e +a e )v q \x z \ 

95 - 82v + I73v 2 - 85v 3 



1 + C F a s 







- 3 ) 


\2v 


(1- 


AtcJ 



101 



(15) 



and 



~AA 



1 + C F a. 



71 

2^ 



7T 3 \ 43 - 30v + 15v 2 + 71v 3 ' 
2 ~ 4tt/ 99^ 



(16) 



Now all ingredients are available to treat the forward-backward asymmetries of heavy 
quarks at QCD one-loop level. The forward backward-asymmetry A FB measures the 
fermion events in the forward and backward hemispheres and is therefore defined as 



da ° da 

dcosO — — / dcosO 



At,, 







dcosO 



dcosO 



( dcos 9 ° + [ dcos 9 
J d cos 9 J 

-1 



d 



(17) 



a 



dcos 9 



which immediately gives with Eq. (Jj]) 



At,, 



Aa v + a L 

For heavy-quark production, the 0(a) corrections for A q FB take a particularly simple form 
when expressed in terms of the effective couplings g % \ with i,j = V,A. Using Eqs. (H) 
and fli~4] ) we obtain 

A q , - fig) 

Close to the Z-peak, the \xz\ 2 propagator terms in the couplings dominate and the cor- 
responding Born expression A 9 FB/0 reduces to the approximate formula given in Ref. JT2 . 



Fig. 2 shows the energy dependence of A q FB for charm, bottom, and top quark 
in the Born approximation (solid line) and including the 0(a) final state corrections 
(dashed line). For the running of the strong coupling we have used the central value 
af\M z ) = 0.118 with five active flavors on the Z-peak. Crossing the top-quark threshold 
the appropriate matching condition has to be imposed. For charm and bottom quark, the 
fixed-point masses m c (m c ) = 1.3 GeV and m^m^) = 4.33 GeV correspond to m c (M z ) = 
0.78 GeV and m^iMz) = 3.20 GeV, respectively. In the top-quark case, we have chosen 
M t = 180 GeV as pole mass, which gives the fixed-point mass m t (m 4 ) = 172.1 GeV. 
In each individual plot of Fig. 2, the dashed-dotted line refers to the right-hand scale 
measuring the absolute difference between Born and 0(a) predictions, i.e. AA q FB = 
A q FB/1 —A q FB/0 . For certain center-of-momentum energies a subtle cancellation of the 0(a) 
expressions in the numerator and denominator of Eq. ( |T%| ) takes place and gives dynamical 
zeros for AA q FB . 

Fig. 2(a) highlights the two specific energy values E cms = 6.58 GeV and E cms = 
89.362 GeV at which the massive one-loop QCD corrections vanish. At the first energy 
point, we find A C FB/Q (6.58 GeV) = A C FB/1 (6.58 GeV) = 0.0038. Of particular interest is 
the energy domain close to 89.362 GeV, where both lowest-order as well as 0(a) one-loop 
contributions to the forward-backward asymmetry vanish. These cancellations originate 
from the dynamical interplay of the involved electroweak couplings at Born level. Note 
that the 0(a) result factorizes with the Born term so that the isolated one- loop zero 
occurs only for the lower energy value E cms = 6.58 GeV. 
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As illustrated in Fig. 2(b), a similar situation arises for the bottom quark. Here, 
our theoretical predictions for the dynamical zeros are A^, (23.32 GeV) = —0.1015 with 
AA* B (23.32 GeV) = 0, and A* B/0 (85.235 GeV) = A* B/1 (85.235 GeV) = 0. 

An interesting observation is that for charm and bottom quark the lower energy 
values yielding dynamical zeros correspond both to approximately the same mass scale 
v = 0.9512. Note that on the Z-peak the predictions of Figs. 2(a) and (b) agree very well 
with the recent estimates of Ref. 



Due to the high production threshold of the top quark, only one first-order a s zero 
is located at 848.80 GeV in Fig. 2(c). The corresponding forward-backward asymmetry 



at Born level is A t FB/Q = 0.5557. 



To isolate the energy dependence of the non-trivial one- loop contents in A 9 FB/1 , we 
introduce the function $ 9 (t> ) in the following manner 



A q - A q 

n ra/l n FB/0 



n/ 

1 + — $ q (v) 



(20) 



Note that Q q (v) contains additional couplings g vv and g AA , which depend solely on v once 
the quark type q is fixed and the corresponding running mass m q (q 2 ) is implemented. 

In the following, we shall use the notation Q q for the exact 0(a) result derived from 
the analytic expressions for a F and a T as given in Eq. (||) and in Ref. jlT], respectively. 
Any additional subscript will indicate approximate representations different from the exact 
formula $ 9 . 

The exact result $ 9 involves all the analytic integral solutions listed in Ref. 
yielding complicated and lengthy expressions. However, using simple polynomial inter- 
polations for the couplings g % \ we can obtain compact and very accurate approximations 
that allow for an easy implementation of the forward-backward asymmetry. For example, 
the substitution of the 0() expressions Eqs. (|T3"D, (|T5|), and fll6|)) into Eq. (^) gives such 
a useful representation, which we denote as $|. 

Further simplification is accomplished by expanding numerator and denominator 
of Eq. (|H3|) up to first oder in a s . After inserting the numerical values for the relevant 
electroweak couplings and the running of the quark masses, we find the following lowest- 
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order (L) interpolations 



^ c (v) = -191.175 + 413.338 w- 223.098w 2 , (21) 
= -0.427 + 5.187w- 5.070 v 2 . (22) 



Note that in this representation $> L is practically insensitive to the distinct couplings and 
masses of bottom and charm quark, i.e. $^' c (l) ~ —1. Only the top quark requires a 
separate parametrization due to its exceptionally high mass. 

To illustrate the high accuracy of these compact formulas, we have plotted in Figs. 3 
the representations $3 and $ i; and the exact 0(a) result $ for charm, bottom, and top 
quark. Furthermore, we have also included the following approximation obtained from an 
expansion of the form factors to second order in y/£ (first order in to) as given in Ref. || 



*l = -! + 1 y/i+k 7+ ±tt 2 - 2 In ( +ln 2 ( ^ 



-H- ' V / J - (23) 

+ % 

Note that $ z contains a part which explicitly depends on the electroweak couplings. This 
formula has primarily been used for approximate estimates on the Z-peak. 

Comparing Figs. 3(a) and (b), we recognize that in the charm- and bottom-quark 
case for v > 0.8 all specific properties as flavor and mass effects fully factorize with the 
Born contribution A FB/0 , and $ exhibits the same universal functional dependence. In 
the region v > 0.9, which corresponds to center-of-momentum energies above 5 and 20 
GeV for charm and bottom quark, respectively, <3> L gives the best representation. Closer 
to threshold, $3 provides more accurate numerical estimates. 

Fig. 3(c) demonstrates that for the top quark $ L and $ 3 yield far better results than 
$ z . Higher-order mass terms in $ z should be taken into account to reach a comparable 
precision. 

To complete this discussion, we present in Tabs. 2-4 the explicit numerical values 
of Born and 0(a) contributions to the forward-backward asymmetry for charm, bottom, 
and top quark with the same choice for ct s (M z ) and the fixed-point quark masses m q (m q ) 
as in Fig. 2. In the fourth column the impact of the 0(a) contributions on the Born 
result is given by (A% B/1 - A q FB/0 )/ A q FB/0 in percent. Again, A q FB/3 , A q FB/L , and A q FB/z 
denote the estimates stemming from $3, and respectively. For these approximate 



results we also provide in the adjacent columns the relative error to A q FB/1 in percent. By 
comparing these relative errors with the corresponding percentages capturing the exact 
0(a) modifications, one obtains a reliable measure of the quality of the approximations. 

In this work, we have investigated the QCD one-loop corrections to the single-jet 
forward-backward asymmetries for massive quarks off the Z-resonance, which modify the 
Born predictions by approximately 3%. After finding compact Schwinger representations 
for the C-odd contribution in the angular distribution of the produced quark, we derived 
simple polynomial expressions that provide excellent approximations for the 0(a) forward- 
backward asymmetry. 

Due to a dynamical cancellation of C-even and C-odd components in the single-jet 
asymmetry, 0(a) zeros occur close to the Z-peak at 89.362 GeV for charm production 
and at 85.235 GeV for bottom production. A thorough experimental investigation of this 
effect should reveal interesting new insights in the physics at higher-order QCD level and 
might even allow for a detection of signals beyond the Standard Model. 
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Figure Captions 



Fig. 1: Schwinger representations for the QCD one-loop corrections to a F in the differ- 
ential rate of the production process e + e~ — > 7, Z — > qq(g). The exact result 
is compared to the polynomial representations of degree m = 2, 3, and 4. 

Fig. 2: The single-jet forward-backward asymmetry A q FB in the Born approximation 
and at 0(a) as a function of the center-of- momentum energy for (a) charm, 
(b) bottom, and (c) top production. The dashed-dotted line refers to the right 
ordinate where AA q FB denotes the difference between Born and 0(a) results. 

Fig. 3: Non-trivial one-loop contents of A q FB/1 as a function of v = ^1 — Am 2 , / q 2 for 
(a) charm, (b) bottom, and (c) top quark. The exact result $ 9 is compared to 
the representations $3, Q L , and $ z . 
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Table Captions 



Tab. 1: Coefficients for the lowest-order Schwinger representations of a F . 

Tab. 2: Forward-backward asymmetry for e + e~ — > 7, Z — > cc at Born level and up to 
0(a) compared to estimates given by <3> 3 , $ L , and $ z (a s (M z ) = 0.118 and 
m c (m c ) = 1.3 GeV). 

Tab. 3: Forward-backward asymmetry for e + e~ — > 7, Z — > 66 at Born level and up to 
0(a) compared to estimates given by $3, $ L , and $ z (m b (m b ) = 4.33 GeV). 

Tab. 4: Forward-backward asymmetry for e + e~ — > 7, Z — > tt at Born level and up to 
0(a) compared to estimates given by $3, $ L , and $ z (M 4 = 180 GeV). 
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Figure 3(a) 
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Figure 3(b) 
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Figure 3(c) 
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Table 3 





At 

FB/O 


4* 

n FB/l 


% 


At 

FB/3 


% 


4* 

FB J L 


% 


At 

FB/Z 


% 


350 
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